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Abstract. The common envelope event is one of the most important and uncertain 
evolutionary stages that lead to formation of compact binaries. While the problem is 
almost 30 years old, its theoretical foundation did not progress much from the first pro- 
posed consideration. For many years, the simple estimate provided by avl-formalism 
has been intensively used by population synthesis studies and, not surprisingly, fre- 
fV| quently contradicted observations. In recent years, the advancements in our studies 

r/y of stellar structure, progress of the numerical techniques for hydrodynamical simula- 

tions as well as increase of the computer power and new observations started to bring 
improvements to our understanding of the common envelope phase. We review main 
physical processes taking place during the common envelope phase from the theoretical 
O ■ point of view and how they affect the values of classical formal parameters. In partic- 

V- 1 ' ular, we discuss the energy budget problem - what are the energy sources, sinks and 

c/3 \ what is the condition for the envelope to disperse, as well as the importance of choosing 

C$ . the definition of the remnant core to the common envelope outcome. 
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1. Introduction 

A key event in the formation an d evolution of most compact bin aries is the short-lived 
common envelope (CE) phase ( Qstrikerl 1 19761 : IPaczynskil fl976 ). The CE phase could 



start when the largest of the two components (the donor star), or both components at the 
same time, had initiated dynamically unstable mass transfer to its companion star (see 
also §2). This occurs, e.g., if the donor's envelope expands faster than its Roche lobe, 
and eventually engulfs the companion star. The drag imposed by the CE transports both 
energy and angular momentum out of the orbit and thus shrinks the binary. 

Depending on whether the envelope can be dispersed sufficiently fast, the system 
will end up either as a close binary or a merged single s tar. The standard dis criminator 
test is known as the energy formalism (Web blnHl98ilLivio & SokeT l988). in which 
the energy difference between the orbital energies before and after the event is compared 
with the energy required to evacuate the envelope to infinity. The latter, denoted £bind 
in the equation below, is commonly taken to be the sum of the potential energy of the 
envelope and its internal energy. 

Gni\ni2 Gm,\^ni2 

-Ebind = A£orb = £orb,i ~ £orb,f = ~ v (1) 

La\ Lai 

Here a\ and af are the initial and final binary separations, m\ and ni2 are the initial star 
masses and m\ c is the final mass of the star that lost its envelope. 

1 



2 



Ivanova 



'bind 



J a 



surface 



(*¥(m) + e(m)) dm 



(2) 



Here *P(m) = -Gm/r is the gravitational potential and e is the specific internal energy. 

For simplicity of use, specifically for population synthesis simulations, a simpli- 
fied parametrization was introduced. The binding energy is expressed then in terms of 
the donor envelope central concentration, using a parameter A: 



'bind 



AR\ 



(3) 



Here m\ fi -m\ — m\ )C is the mass of the removed giant envelope and R\ is the radius of 
the giant star at the onset of CE. 

To measure the energy transfer efficiency from the orbital energy into the envelope 
expansion, another parameter, a, was introduced. The balance of energy is written then 
as 



a A 



Gm\ iC ni2 Gm\m2\ Gm\m\ 
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(4) 



Since for many stars (at least low-mass star, as detailed models show) A = 1 and a < 
1 due to energy conservation, it is very common, if detailed stellar models are not 
available, to accept a A = 1 . 

Another a pproach is to consider, instead of the e nergy balance, the angular mo- 
mentum (Neleman s et al.ll2000l : iNele mans & Tout 2005). It is known as y-formalism: 



A/lost Ji - Jf 



Ji 



y 



mi, 



(5) 



Here 7; and Jf are the angular momenta of the initial and the final binaries. 

y-formalism could be used to explain the formation of the systems for which it 
was originally calibrated with observations (double white dwarf binaries). However, 
it has been shown that this reverse engineering produces a range of values of the pa- 
rameter y whic h ca n not be extended to the case of the CE in an arbitrary binary (see 
Webb inM 120081 and lWoods et all I2010L for more formal mathematical explanation see 
Woods et al.1 1201 lal) . It has to be appreciated as well that even in the case of dou- 
ble white dwarf binaries y-formalism implies an energy generation during the first 
episode of the mass transfer, reconstructing as suc h a non-dynamical e vent but not a 
CE (likely, a stable mass transfer, as was shown in IWoods et al.ll201 lbl) . More recent 
comparisons of the observed populations of post-CE binaries with t he population syn- 
thesis models had also shown severe inconsistency with y-formalism ( Davis et"aT1l2010l : 
Zorotovic et al.ll2010h . Considering recent theoretical and observational disagreements 
with the y-formalism statements, this formalism, in its current form, can not be recom- 
mended for studies of arbitrary CE events. 

Let us review in more detail the energy formalism. There, the simple parametriza- 
tion with aA = 1 had also started to fail once more statistics from observations came. 
In particular, to explain the observed f ormation rate of low-mass X-ray binaries, 1 per 
mln yr per Milky Way (Romani 19981) . it is apparently necessary to adopt unrealisti- 
cally high values of a (IKalogeralll999l) : e.g., lYungelson et al.l d2006l) showed that the 
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Figure 1. Main phases of the common envelope. Showed the case of a 1.6 M Q 
red giant and a 0.3M Q WD. Dashed line s show the lines of constant mass. Used data 
from Id hydrodynamical simulations in llvanoval d2002l) 



observations could be matched with a\ > 2. The situation became even more dra- 
matic with the use of values of A taken from detailed stellar models: in massive giants 



A <sc 0.1 (ITauris & Dewi 2001; Podsiadlows ki et al.l l2003). It has been shown as well 



that with a < 1 only an intermediate-mass companion could avoid a merger; this chal- 
lenges the formation of a low-mass X-ray binary with a low-mass companion in general 
(lJustham et al.ll2006h . 

The situation may seem dire, but progress in several subproblems has been made. 
In order to understand the roots and appreciate the profundity of the problem once the 
straw-clutching of parametrization - which is done by forcefully fitting the population 
synthesis runs with observations - are discarded, one has to resolve two seemingly 
obvious questions: 



What is really A? 
What is really a? 



2. Phases of the common envelope event 

For further understanding, it is important to realise that a CE event is not a single self- 
similar and purely dynamical event for all the binaries, but most likely a seq uence of 
several phases, where each operates on its own timescale (IPodsiadlowskill200l[ see also 
Fig. 1): 

• Loss of corotation 



The start of the spiral-in could be caused by, e.g. 
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(a) a dynamically unstable (runaway ) mass transfer (for a linear stability analy- 



sis of the mass tran sfer see, e.g., Hiellmi ng & Webbink!ll987l : ISoberman et al 
ll997l : lRitteri [T999) 



(b) 



(c) 



a dynamical inst ability (the Darwin's instability, Darwinlll879 ), or a secular 
tidal instability (fHu3ll98dl : ILai et all 1 19931 : lEggleton & Kiseleva-Eggleton 
l200lh . when the spin angular momentum of the system is more than a third 
of its orbital angular momentum 

in some special cases, e.g. in nova systems because of the expansion of the 
nova shell engulfing the companion, or in systems with a compact compan- 
ion where a stable but very fast mass transfer leads to a formation of the 
extended envelope around a comp act object: e.g., reincarna tion of a white 
dwarf in the form of a red giant dNomoto et al.lll979l l2007h . or "trapping" 
the accreting material in the form of the envelope around a neutron star in 
case of supereddington mass tr ansfer dBegelmanll 19791 : iHouck &~C hevalier 
199ll : lKing& Begelman| [T999h 



The loss of the corotation occurs on a dynamical timescale, although the stel- 
lar structure at that moment is strongly affected by mass-transfer history before 
the dynamical instability sets on, and this stage could last hundreds of years 
dPodsiadlowski et al.ll2002h . 



Plunge-in 

A rapid spiral-in, during which all the frictional energy is deposited in the enve- 
lope, drives its expansion and may lead to its dynamical ejection. This stage is 
purely dynamical. 

Slow spiral-in 

When the orbital decay transferred enough energy to cause significant envelope 
expansion, the spiral-in slows down. This is a self-regulating phase when all re- 
leased frictional luminosity is transp orted to the surface where it is radiated away 
(Meyer & Meyer-Hofmeister 1979). This phase is non-dynamical and operates 
on the thermal time-scale of the envelope. 

Termination 

Slow spiral-in ends with either the ejection of t he envelope (e.g., via delayed 
dynamical ejection, Iivanovall2002l lHan et al.ll2002h . or by a (slow) merger. 



Further evolution of the survived binary could also be affected by a circumbinary 
disk, winds from the giant's remnant and so on. 

The dynamical plunge-in phase, when considered as a separate "e ncapsulated" 

event , is the best studied pha se of the CE to date (for a thorough review see iTaam & Sandquist 
20001 : Taam & Ricked [2010L and references therein). This stage proceeds dynamically 
and no additional energy sources or sinks affect its evolution. The com plications arise 
from: (i) this stage usually does not end with a dynamical ejection dRicker & Taaml 
20081 : ITaam & Rickeill20 10). but instead hydrodynamical simulations show the transi- 
tion into a slow spiral-in; (ii) the initial conditions are not properly modeled and this 
may affect the final results (see more details in §5). 
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3. Story of the /1-parameter 

A is the factor that links the "properly" calculated binding energy from detailed stellar 
models and its parametrized form. To be exact, by definition it is: 



Gm x m u f surface 



f 



QV{m) + e(m)) dm = Q sm + U em = W env (6) 



Here f2 env is the potential energy, U em the internal energy and W env is the total energy 
of the envelope. If a code that calculates detailed stellar models is available, and if 
the core-envelope boundary is known exactly, A can be calculated for any star at any 
evolutionary stage. 

Intuitively, A was meant to be a measure of central concentration of the star, which 
would translate into the expense side of the energy budget. 

3.1. Simplification of detailed stellar models 

Is some works, an over-simplification to obtain A from detailed stellar models is made, 
by connecting the potential energy and the internal energy as 

— ^bind — W^env — ^env + U em/ — 2^ env 

This, presumably, implies the use of the virial theorem. It should be appreciated that 
this expression is not valid for the following reasons: 

• Virial theorem indeed states that, for a star in hydrostatic equilibrium, 2K+Q, = 
(here K is the total kinetic energy), however this is valid only for a the entire 
star. 

• K = U only in the very specific case of constant adiabatic exponent Ft, = 

(fllip) d = ^ n ^ s * s not at a ^ not tne case * n e, & - ' P^^l ionization lay- 

ers, or in massive stars. E.g., in massive stars, Ft, — > 4/3 and, since K — > U/2, 
then £bind = W — > 0, which is far from ^bind = -^env/2 as in the Eq. $7}. 

If one defines A g as only dependent on the potential energy (as in iTauris & Dewil 

1200 lh . and so A s /A = £'bind/^env> then simplification as in the Eq. (0 would predict 
Ag/A = 1/2; however detailed stellar models, as expected, show a very different ratio 
(see example on Fig. 2). 

3.2. A constant value: how good is that? 

In the case when detailed stellar models are not available, the most commonly accepted 
value is A = 1. This value is rather good approximation for low-mass red giants, espe- 
cially with masses ~ 2M Q . In a general case, however, A is the functi on of mass and 



radiu s, and it changes as the star evolves. For mo st stars. 0.1 < A < 45 (IDewi & Tauris 
120001) . and massive stars could also have A < 0. 1 (IPodsiadlowski et al.l 120031) . 



The Eq.© tells that the adoption of A = 1 leads to about the same error in the final 
binary separation by the same factor as a "properly" found A differs from 1. Tabulated 
values o f A for stars as a function o f mass and radius, and /or fi tting formulae, can be 
found in lvan der Sluvs et all (l2010h . lLoveridge et all feOlCh and lXu & Lil fcOlCh . 
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Figure 2. i as a function of mass shown on example of 20 M star when it has 
R = 950/? Q (Z = 0.02, overshooting 0.2 of the pressure scale and no wind loss). For 
comparison A as in Eq.(2) (thick solid line) and A g (thin solid line, see discussion 
§3.1) are shown. Dotted lines correspond to several possible core definitions, as 
discussed in §3.3, 3.4. 



3.3. Uncertainty given by the boundary problem 

Even with detailed stellar models, we still need to decide on where to p lace the bound- 
ary be tween the core and the envelope. It has been first noticed in iTauris & Dewi 



(2001), that placing this boundary within the hydrogen shell between the He core and 
the bottom of the convective envelope, A can change by a factor of 10-70, depending 
on the star, where the discrepancy is more marked for more massive stars. This is ex- 
plained by the fact that the binding energy of the H shell far exceeds the binding energy 
of the convective envelope above. 

The uncertainty in the value of A due to an ambiguous core definition, as discussed 
in §3.2, translates directly to a similar uncertainty in the final binary separation. It is 
clear therefore that the ambiguity in core definition leads to the same or even larger 
uncertainty than the use of a constant A for all the stars. Worse still, this uncertainty is 
invariably overlooked in the studies of CE outcomes, and the accepted definition of the 
core in a particular set of simulations is usually the least clearly discussed condition. 
Note that the vagueness in placing the core boundary also comes from the fact that the 
remnant of a giant in the formed binary (immediately after the envelope ejection) does 
not have to be the same as the post-CE remnant we observe today, as the remains of 
the Hydrogen shell could be later removed by strong winds, similar e.g. to Wolf-Rayet 
winds. 

A variety of possible definitions of where is the boundary between the remaining 
core and the ejected envelope could be found in literature: 



• the minimum possible core mass is at X = (bottom of H shell) 
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at the maximum nuclear energy generation within H shell dTauris & Dewill200lT) 



at the maximum nuclear energy generation plus the condition on the thickness of 
the remaining envelope, which itself is a function of the evolutio nary status of the 



donor (for low mass red giants and asymptotic giant branch stars. iDe Marco et al. 
20lTh 



wher e the nuclear energy generation falls below some threshold (IDe Marco et al. 



2011) 



the central mass which contains less than 10% hydrogen ( Dewi & Taurisl l2b00) 

core is everything below the location where X = 15% ( Xu & Lill2010h 

where d 2 log p/d 2 m = within H shell (IBisscherouxll 19981) 

where the function of the binding energy y = sinh'^isbind) has the tran sition be- 
tween a sharp increase and a fairly slow increase in the outer envelope dHan et al.l 



1994) 



where the entropy profile has a transition between its increase and flat part dTauris & Dewi 



2001) 



where the effective poly tropic index has discontinuity (Hjellming & Webbink 



1987) 



the bottom of the convective envelope (top of the H shell) gives the maximum 
possible core mass, as convective envelopes are known to expand upon mass 
loss. 



Last three definitions describe the same place, the transition between the radiative 
zone of H shell and the convective envelope. Note that not all definitions are applicable 
to every star. Some are created only to low-mass giants or asymptotic giant branch stars 
(e.g., thickness of the remaining envelope), some conditions can not be either found or 
determined uniquely in all the stars (e.g., the condition on d 2 logp/d 2 m = can not be 
determined in case of 20M G as on Fig. 2). Unfortunately, none of the above definitions 
has a solid physical foundation explaining the choice. 



3.4. Core boundary consideration via post-thermal readjustment 

Anoth er method to define the core remaining after a CE event was proposed in llvanoval 
(120111) . The idea is that different parts of the H shell react differently to the very rapid 
mass loss; both immediately after the envelope ejection and subsequently on the thermal 
timescale. Indeed, the bifurcation in the reaction of H shell material - whether the shell 
re-exp ands after the mass loss or collap ses - has been fo und long ago for low-mass 
giants dDeinzer & von Sengbuschlll970h . In llvanoval ( 2011 ) it was shown that similar 
divergence points exist in the H shell of more massive giants as well. If mass is removed 
below such a divergence point, the remnant contracts on thermal timescale from its 
post-CE state, if the mass is left above such divergence point, the star expands during 
its thermal readjustment; it could experience strong thermal pulses, or even develop 
convective envelope. The subtle difference with the low-mass case is that, immediately 
after the fast mass removal, the remnant, as it contains a non-degenerate core, always 
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starts to expand compared to its pre-CE state (similar e xpansion for donors w ith non- 
degenerate cores was found, e.g., in adiabatic studies by lDeloye & Taam ll2bl0). 

The divergence point above coincides with the location where P/p reaches its max- 
imum value within the H shell, so called "compression point" m cp . In general, this 
divergence point does not coincide with any of the ones listed in §3.3. 

Finding the core mass via m cp provides an upper bound on the maximum post- 
thermal readjustment core mass as it is not possible to completely rule out that the 
dynamical phase does not remove mass below m cp . Such extra mass removal however is 
not likely to happen if the characteristic orbital evolution time during the final stages of 
the spiral-in is comparable to the core thermal response time near m cp point - which is 
as short as 10-100 years. If the event proceeds on this timescale, then the core material 
below m cp starts to contract and precludes further mass removal. 

On the other hand, all the mass that could have been left above m cp is guaranteed 
to be pushed out by the thermal timescale mass transfer during star's readjustment. 
If the envelope was ejected during a dynamical plunge-in phase and the core mass is 
nipost-d > m cp , then the core is reduced to m cp without any additional external energy 
source like orbital energy, but only using core's own thermal energy. 

The above property of m cp can be stated as an alternative definition for the bifur- 
cation point: the energy expense required to shed the envelope down to m cp is minimal. 
It is fully reasonable to remove less of the envelope (and have a bigger post-CE mass) 
once £bind < aEorb- The thermal readjustment phase will then remove mass down to 
~ m cp . However, energy-wise, it is not advantageous to remove the envelope deeper 
than ~ m cp . 

The complete understanding what exactly is A (and in other words what is the energy 
required to remove the envelope) however is linked with how the envelope ejection 
happens, on which timescale and what energy can be used during the CE, and therefore 
is amalgamated into the ^-parameter problem. 



4. Story of the a-parameter 

a is the parameter that was introduced to relate the available energy source (orbital 
energy) to the energy required to disperse the envelope, where the latter is assumed 
without questioning to be the binding energy. The energy balance is then written as: 

aAE olh - £ bind (8) 
Energy conservation restricts < a < 1. 

4.1. What energy is really required? 

The original assumption that the energy requirement equals .Ebind carries one of the two 
tacit assumptions: that an envelope is dispersed once its W em > 0, or that an envelope 
with W em > is unstable. This approach has been questioned in ilvanova & Ch aichenets 



(120111) . for the following reasons: 



• a star with W > can be "kinetically stable" (i.e. additional energy is re- 
quired to overcome an energy barrier between the bound and unbound states, 
Bisnovatvi-Kogan & Zel'DovicrJ[l967h . 



Common envelope: the progress and the pitfalls 



9 



a star or its part (e.g. envelope) is unstable against adiabatic perturbations if the 



first a diabatic exponent Ti is approximately a constant throughout the star (IChiu 
1 1968b and Ti < 4/3. This condition however is not the same as W > 0, which 
requires a constant third adiabatic exponent < 4/3 (Hans enet all 12004]) . The 



relation between the adiabatic exponents is Ti = 1 + TiV a d, where V a( j is the 
adiabatic temperature gradient. 

The energy conservation equation in La grangian coordinates applied to each mass 
shell in the envelope is written as (see, e.g.. IColgate & White1ll966l : iKutter &"S parks 



1972): 



d (u 2 \ d(l/p) 

Here u is the velocity, P is the pressure and p is the density. When heat Q (e.g., fric- 
tional energy from orbit's dissipation) is imposed into the envelope with some arbitrary 
distribution 8q{m) > (Sq(m) is per mass unit), such that 



J cc 



surface 

8q(m)dm = Q, (10) 



the energy conservation for any given lagrangian shell can be rewritten as a version of 
the Bernoulli equation: 



(6q(m) + x ¥ + e + -) sta rt = (-" + ¥ + e + -)exp = 2 = const (11) 



PI P 

-)start = (-U 2 + ¥ + € + -] 
p 2 p 

A Lagrangi an shell then will reach the point of no return in its expansio n when its 
I > 



igrangian shet l men will rea ch the point oi no return in its e xpansion when its 
(e.g.. iBisnovatyi-Kogan & ZerDovichlll967l : ISparks & Kutterill972h . Past sta- 



bil ity analyses showed that a star starts outflowing once its envelope obtained positive 
S ( Bisnovatyi-Kogarj|2002l and the references therein). Same will occur in the case a 



part of the common envelope has obtained positive E, and the whole common envelope 
should fully disperse if the entire envelope would have positive X: 



/-•surface / P(rn)\ 

Q + rP(m) + e(m) + — — \dm > . (12) 

Jcore \ P(m) ) 

The development of such quasi-steady surface outflows operates on the timescale 
on which the envelope redistributes the dumped heat, i.e. the thermal timescale of the 
envelope, or the same timescale on which the self-regulating slow spiral-in operates. 
These outflows are not likely to take place during a dynamical plunge-in phase. It is 
important, that due to the presence of a non-negative term P/p, this condition occurs 
before the envelope's total energy become positive and as such eases the survival of 
a binary with a smaller by mass companions. In particular, the use of this criteria in- 
stead of the original energy for malism allows the formation of low-mass X-ray binaries 
(llvanova & Chaichenetsll201lh . 

4.2. Other energy source 

In studies of binaries that, by the energy formalism, are doomed to merg e, it was found 



that the nuclear energy source c ould play a role in the envelope ejection (llvanoval 2002; 



Ivanova & Pod siadlowskil l2003l) . 
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Such binaries enter the common envelope as any other binary, but do not merge 
during their dynamical plunge-in phase, proceeding to the self -regulating spiral-in phase. 
If the companion is a main sequence star, then at the end of the slow spiral-in it starts to 
overfill its Roche lobe while the core of the giant is still well within its Roche lobe. The 
mass transfer is forced by continuing orbital shrinkage due to drag from the envelope, 
its rate is increasing and, depending on secondary and the mass ratio, is from 0.01 to 
100 M Q per yr. A thick stream of hydrogen-rich m aterial is able to pen etrate deep to- 
wards the core, reaching even the He burning shell dlvanova et alj |2002). The outcome 
of the steam penetration depends on the initial stream entropy. The latter depends on 
the mass of the secondary: a lower mass main sequence star has a lower entropy, and 
its material is able to be injected deeper in the giant core. 

In some cases, the outcome of such rather slowly proceeding merger is the com- 
plete explosion of the He shell. The released nuclear energy during explosive hydrogen 
burning could exceed the binding energy of the He shell, which is in massive stars can 
be about few 10 51 ergs. The rest of the common envelope has much smaller binding en- 
ergy and is ejected during the same explosion as well, leaving behind a compact binary 
consisting of the core of the giant and whatever remains from a low-mass companion 
after the mass transfer. The effective a, compared to the available orbital energy, is up 
to 100. 



The main consequences of this explosive CE ejection are (IPodsiadlowski et al. 
20101) : 



• a less massive companions could survive the CE, this helps the formation of low- 
mass X-ray binaries and could explain their observed formation rates 

• a fast rotating, hydrogen exhausted core is formed. This core could be a long y- 
ray burst progenitor and Ic type SN progenitor at the same time, explaining their 
observational connection. 



4.3. "Useful" energy 



One of the important questions is which fraction of the available internal energy can 
be used to eject the envelope. Arguments both for dHan et al.lll994l 1 19951 . 12002D and 
against dSokerj|2002h ISoker & Harpazll2003r) using the recombination energy have been 
advanced. The question is not yet settled, ev en though it has very important conse- 
quences discussed in details in dWebbinkll2008r ) : if this energy should be included, then 
some of common envelope events begin with positive total energy of the giant envelope, 
bringing the effective value of A -parameter below zero ! 

We note in addition that recom bination energy leads to yet another attempts of 
parametrization in some works: e.g. lHan et al] (1995) introdu ces an "efficienc y of use 
of internal energy" and thus special second a-parameter, a^. Web binkl d2008l) instead 
proposed to define an ejection efficiency /? that has the desired property, < /? < 1: 



P 



A£o rb - U e 



(13) 



The latter definition however only serves to hide the unsightly blemish of a not falling 
within the range < a < 1. Overall, attempts to parametrize the situation with some 
sort of a new parameter instead of simply writing the full energy balance could only 
confuse the reader. 
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4.4. Energy losses 

The energy balance between the available and required energies as in the Eq.(l) implies 
that all the heat is transferred into the expanding envelope and that dispersed envelope 
has zero velocity at infinity. This is however far from the results obtained in detailed 
simulations: 



Radiative losses: 

If the CE did not end during the dynamical plunge-in phase, it proceeds with the slow 
self-regulating spiral-in. In this case, most of the orbital energy dissipates at those 
shorter orbits, when the envelope is already expanded (see also Fig. 1). While the giant 
donor star has expanded, its evolves along the giant track toward higher luminosities 
and keeping almost the same effective temperature. Its luminosity L is then increased 
as oc R 2 , while the nuclear energy sources do not catch up with the increased L. It is 
the orbital energy that is converted into internal energy and then is radiated away, as 
for all stars evolving on thermal timescale. It is possible to loose this way a significant 
fraction of the available orbital energy. E.g., we find that a binary consisting of 20 
and 5 M Q companions is loosing during the slow spiral-in 40% of the released orbital 
energy by radiating it away from the surface of the CE, in case of 1M Q companion the 
slow spiral-in stage lasts very shortly and not much energy is lost, but in the case of a 
1.6 M Q giant and a 0.3 M Q companion the losses are close to the 100% of the released 
orbital ener gy until the del ayed dynamical ejection starts (data is from the simulations 
described in ilvanovall2002h . 



Non-zero velocity at infinity: 

It is currently adopted for simple analysis that at the infinity velocity of the expanding 
envelope u - 0: this is how the full energy conservation equation Eq. (fTTb is translated 
into simplified energy formalism like in the Eq. (fT2l or E. (1). However, hydrodynam- 
ical simulations performed with use of a smooth particle hydrodynamics, show that 
the kinetic energy of the ejected envelope at the infinity could be from 15% to 70% 
of Iffn rh f | (data is from simulations performed for Lombardi et al. 2006; llvanova et al. 
201Ch . Values of the kinetic energy do differ depending on the mass ratio of the inter- 
acting stars and also on the core to mass ratio in the interacting giant. 

The consideration of just two possible ways to loose energy compared to the pure con- 
version of the available energy into only the envelope ejection energy shows that the 
effective a (meaning that all possible energy sources are also taken into account, not 
only the orbital energy) is not only less than 1, but also is a function of the stellar 
masses and mass ratios in the binary system. There is no reason whatsoever to expect 
one value of a to apply to all CE events. 



5. Start of the CE and its influence on a and A 

A not particularly well studied CE phase is the onset on the CE (see Fig. 1). This is 
not an immediate process and involves the evolution of the donor before the dynamical 
evolution stars. The whole list of possible unsolved problems during the initial stage 
when the binary is loosing its corotation includes such questions as the possibility for an 
enhanced wind losses just prior to the Roche lobe overflow, dependence of the critical 
mass ratio, as a condition to start runaway mass transfer, from the giant structure. Even 
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without consideration of this, we know that hydrodynamical simulations experience 
the problem to start the CE. One of the important factors is how far the donor is out of 
corotation with the binary. 

In order to have successful start of the CE, the initial conditions for common- 
envelope hydrodynamical simulations are such that the dynamical plunge-in phase is 
literally forced to start quickly. It can be done by either con sidering a donor being, e.g., 
not in a corotation with a binary (as in iRicker & Taamll2008l) . or by having some degree 
of the Roche lobe overflow initially, or by introducing initial tidal perturbations (some 
variation of non-corotation). Choosing the degree of non-corotation leads to different 
final results for the ejected masses and the masses of formed circumbinary envelopes. 
A fully self-consistent solution in term of initial condition so far is obtained only in case 
of physical collisions between a giant and an intruder , where initial non-corotation is a 
physical sound condition (e.g.. lLombardi et all 1200 6). However these simulations are 
always accompanied by a significant kinetic energy of the ejected envelope at infinity 
and it is not clear that they can be extrapolated to the case of a normal CE. Nonetheless, 
it is clear that the choice of the starting condition for the dynamical plunge-in phase 
affects the final a by taking away different amount of energy in the form of the kinetic 
energy at the infinity. 

Another reason for either A, or the combination of a A to be affected by this not 
well studied phase is that the dynamically unstable mass loss/mass transfer in fact may 
start not immediately. It could have been pre ceded by either stronger than usual winds 
(just prior the full Roche lobe overflow, e.g. lEggletonll2002l) . or the mass transfer be- 
come dynamically unstable only after a significant fract ion of the envelope has been 
transferred on thermal timescale dWoods & Ivanovall201lh . or, if mass transfer is super- 
Eddington, the accretion energy released i n the rapid mass -transfer phase unbinds a 
significant fraction of the giant's envelope dBeer et all 120071) . Effectively, in all three 
cases, the envelope mass is decreased. This reduces how much energy has to be spend 
then later on to expel the formed common envelope. It is not entirely clear with which 
of the two parameters this decrease in mass should be associated. 



6. Conclusions 

As of today, we can list the main pitfalls in studies of the CE using parametrization. 
Those pitfalls should be either avoided, or carefully discussed in terms of the (large) 
uncertainties they provide for the population of post-CE binaries: 

• in general, no parametrization can be fully trusted 

• the parametrization in the form of y-formalism can be used only to replicate 
double white dwarf binaries formation, similar to those used for calibration 

• the parametrization for A can not use a constant value: A always depends on mass, 
radius and evolutionary state of the star, and A value can differ from 1 by up two 
order of the magnitude. 

• value of A depends as well on the choice of the remnant core mass, and with 
different core masses your mileage may vary as well by up to two orders of the 
magnitude 

• Energy sources such as recombination energy can lead to a negative value of A 
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• a deposition into the common envelope of the energy equal to £"bind does not 
guarantee the CE ejection! 

• an assumption < a < 1 is not valid in general case: a can be negative if 
recombination energy should be included, can be a factor of few if the envelope 
is removed by adiabatic outflows, and also can be by up to 100 if nuclear driven 
ejection occurred. 

As we shun the parametrization approach and deal with understanding the under- 
lying physics of the CE event, a much clearer picture emerges: 

• it is likely that we understand now what exactly is the post-CE core, however, its 
mass is not parametrizable in an easy way 

• there is progress in understanding of the ejection condition during the slow spiral 
in. A self-regulated spiral-in results however not only in another criterion to eject 
the envelope, but also in more energy losses, giving as a result a that can vary 
from 0.4 to 3 

• certain binaries could have a very large effective value of a, up to 100, if they 
e.g. experienced an explosive common envelope event 

• recent hydro simulations showed that we are almost able to model a complete 
common envelope event if that includes only dynamical plunge-in phase, though 
a proper onset on the CE (which is an initial condition for the dynamical phase), 
as well as the modeling of self-regulating spiral-in stage is still a challenge 

I cannot stress enough that there is no reason to expect either A or a to have unique 
values for all the binaries, as they do depend on the evolutionary status of the donor, 
the mass ratio, the initial orbital separation - in short, on every parameter and physical 
process that is involved in setting up a CE. Searches done by matching the observed 
post-CE binaries and the modeled population in order to find some unique constant 
value of a might be encouraged only when they are related to well partitioned sub- 
classes of post-CE binaries. A easily varies by two orders of magnitude depending on 
the evolutionary status of the donor and on the adopted core definition, and in some 
cases can be even negative, and available energy ("effective" a) can for some binaries 
exceed the orbital energy significantly. It is unlikely therefore that any value of a con- 
strained from the observations for a group of post-CE binaries can either completely 
rule out proposed formation channels or be extrapolated to other subclasses of post-CE 
binaries. It is certain that further understanding is required for such CE phases as the 
onset on CE and the slow spiral-in stage. 

Acknowledgments. N. Ivanova acknowledges support from NSERC Discovery 
and Canada Research Chairs programs. 
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